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Abstract 

We investigate the unfolding of single polyelectrolyte (PE) chains collapsed by trivalent salt under 
the action of alternating-current (AC) electric fields through computer simulations and theoretical 
scaling. The results show that a collapsed chain can be unfolded by an AC field when the field 
strength exceeds the direct-current (DC) threshold and the frequency is below a critical value, 
corresponding to the inverse charge relaxation/dissociation time of condensed trivalent counterions 
at the interface of the collapsed electrolyte. This relaxation time is also shown to be identical to 
the DC chain fluctuation time, suggesting that the dissociation of condensed polyvalent counterion 
on the collapsed PE interface controls the polyelectrolyte dipole formation and unfolding dynamics 
under an AC electric field . 



1 



I. INTRODUCTION 



To develop new ability to separate and stretch biomolecules, such as nucleic acids and 
proteins, is very important in molecular biology and medical research for its applications 
in single molecule analysis and DNA sequencing [1-4]. Conventional way to determine the 
genetic information of a DNA molecule relies on a series of cumbersome operations of sub- 
cloning and electrophoretic separation [5j. This current method reads the sequence of DNA 
fragments of hundred base-pairs per time, which has become inadequate for individualized 
genetic studies j^j]. To improve the efficiency of sequencing, single- molecule sequencing meth- 
ods have been developed recently (2 - 4 ] , in which the genetic information of individual DNA 
molecules is scanned in a linear fashion. One of the key problems in these methods con- 
cerns the unfolding of DNA molecules to increase the spatial resolution of scanning, as DNA 
molecules are long polyelectrolyte (PE) chains that usually adopt the coiled conformation 
in a typical buffer. In the past decades, many techniques have been developed to stretch P 
chains, including the utilization of optical tweezers [6], fluid flows electric fields 



and other methods 



13- 



12|, 



151 ] . Among these methods, stretching PE chains in electric fields has 



received considerable attention, as electrokinetics has become a viable microfluidic method 
for on-chip molecular analysis jltl \v\ . 

DNA electrophoresis is usually performed in a sieving medium, such as a gel or a poly- 



mer solution 
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191 ] . In that case, DNA chains transform cyclically between an extended 



structure, with one end hooked onto obstacles in the medium, and a coiled structure when 
released from the obstacles 201 ] . If an alternating-current (AC) electric field of appropriate 

chain in a sieving medium can be main- 



frequency is applied instead, the conformation o 
tained constantly in an elongated structure 
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ll| . In a bulk buffer, coiled DNA chains have 



been demonstrated to be stretched in a strong direct-current (DC) electric field, through 
the mechanism of polarization s), 2l|. Nevertheless, there is a drawback to this method: 

22L To 



a strong DC field leads easily to electrolytic dissociation of water molecules 
avoid this drawback, a short-duration pulse or an AC electric field has been suggested, with 
the period shorter than the dissociation reaction time. Moreover, the mean displacement 
of PE chains in an AC field is zero. This feature could be used to trap or to manipulate 
single DNA molecules at some place to facilitate the detection of the sequence, particularly 
when it is integrated with a dielectrophoresis trapping mechanism 17]. Experiments have 
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shown that AC electric field can indeed stretch and orient DNA chains in aqueous solutions 



at low AC 
elevated 
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requency and the chain polarization is largely reduced while the frequency is 



24j . Recently, Wang et al. demonstrated that the conformational transition of 



single PEs possesses a hysteretic nature upon sweeping the frequency of AC electric field [25 ] . 

Without a solid medium to break the symmetry of the coil-stretch transitions, a low- 
frequency AC field should stretch the DNA repeatedly during every half-cycle, in a manner 
of a quasi-DC field. However, unlike DC stretching, the polarization and stretching dynamics 
now becomes important, as the molecule must be fully stretched within one half period. It 
is known that multivalent cations have dramatic effects on the properties and structures 



of DNA molecules in aqueous solutions 26) . Coiled DNA chains can collapse into glo 
particles while the chain charge is neutralized by the condensed multivalent cations [26|, 



DU 



27]. 



Experiments have shown that a collapsed DNA chain can be unfolded in DC electric fields 



if the field strength surpasses a threshold 28J. This threshold E* depends on the chain 
length N, and the dependence has been studied recently by computer simulations 29^32] . 
The results suggested a power law-like relation E* ~ N~ e and the electrophoretic mobility 
of chain depends strongly on whether the chain is unfolded or not. These studies provide a 
solid foundation to separate PE chains by length in free solutions, based upon the drastic 
mobility change when a chain is unfolded by an applied electric field jsi|. Moreover, recent 
study in protein mass spectrometry showed that AC electrospray ionization (ESI) yields 
an effective ionization without fragmenting molecules, with higher ionization efficiency than 
DC ESI 33]. The ionization efficiency is a sensitive function of sample pH and molecular 
conformation. It is still unclear whether the higher efficiency is due to the increased local 
pH or the conformational change or both of the two effects. 

This present work constitutes a fundamental molecular dynamics study of how an AC 
field affects the conformation of a PE chain and the condensation of counterions, and how 
the field strength to unfold the chain varies with frequency, in order to give deep insight of 
the complicated behavior of PE solutions in AC fields. The rest of the paper is organized as 
follows. We describe the model and simulation setup in Sec. II. The results and discussions 
are presented in Sec. III. The topics include: chain unfolding in AC electric fields (Sec. III. A), 
simulation snapshots (Sec. III.B), theoretical explanation (Sec. III.C), time evolution of 
dipole moment and chain unfolding (Sec. III.D), and the study of polarization time and 
fluctuation time (Sec. III.E). A summary is given in Sec. IV. 



3 



II. SIMULATION MODEL AND SETUP 



Our system contains a single linear chain, modeled by a bead-spring chain model. The 
chain is composed of N monomers. Each monomer carries a negative unit charge — e and 
dissociates one monovalent counterion into the solution. The bond connection on the chain 
is modeled by the finitely extensible nonlinear elastic potential 

1 / b 2 \ 

£/fene(&) = -^kb 2 max In M - — \ (1) 

where b is the bond length between two neighboring monomers, b max the maximum exten- 
sion, and k the spring constant. The condensation of the PE chain is induced by adding 
(3:l)-salt into the system. The salt molecules dissociate into trivalent cations (also called 
"counterions" ) and monovalent anions (also called "coions") in the solution. All these par- 
ticles, including monomers, counterions, and coions, are modeled as Lennard- Jones (LJ) 
spheres, described by the potential 



Ae 



LJ 



(,/r) 12 - (<r/r) 6 + 1/4 ,forr<^2a 
^Lj(r) = { L J (2) 

, for r > 



where a and £lj represent the diameter and the hardness of the LJ sphere, respectively, 
and r is the distance between two particles. Since the interaction between the monomers is 
purely repulsive, this setup corresponds to a good solvent condition. Particles also interact 
with each other via Coulomb interaction, which read as 

U c (r) = k B T.^R (3 ) 
r 

where k B T is the thermal energy and Zi is the charge valence of the ith particle. The 
Bjerrum length A# = e 2 /(47reoe r /cBT) is the distance between two unit charges, at which 
the electrostatic energy is equal to the thermal energy, where eo is the vacuum permittivity 
and e r is the relative dielectric constant of solvent. The solvent molecules are not modeled 
explicitly in the study. However, their effects are incorporated implicitly by the following 
three ways: (1) dielectric screening of charge by e r , (2) friction force, —rriiQvi, due to particle 
moving through the solvent, (3) stochastic force, fji(t), owing to thermal collisions of solvent 
molecules. The equation of motion, also known as the Langevin equation, thus reads as 

rriifi = -midfi + F c + ZieE(t)x + ffi (4) 



where m\ is the particle mass, is the friction coefficient, and F c = —dU/dfi is the 
conservative force. In this equation, the temperature T is controlled by the fluctuation- 
dissipation theorem: (fji(t) ■ ffj(t')) = 6k B Tmi(iSijS(t—t'). The system is placed in a periodic 
rectangular box. The AC field E(t) is a square wave applied in x-direction. The period of 
the square wave is ^ and the field strength is \E\. Particle-particle particle- mesh Ewald 
sum [34] is used to calculate Coulomb interaction. 

We simulate aqueous PE solutions at room temperature (T = 300K) and assume that all 
the particles have the same mass m and LJ parameters a and We set Slj = 0.8333fceT, 



k = 5.8333k B T/cr 2 , b max = 2a, \ B = 3cr, and Q = l.Or -1 where r = aym/ '{k B T) is the time 
unit. Since \ B is 7.1 A for water, the length unit a corresponds to 2.4 A. The time unit r 
is about 1.5 ps if m is taken to a typical monomer mass of the order 100 g/mol. This set of 
parameters can be used to model flexible PE chains in water, such as sodium polystyrene 
sulfonate. The dimension of our simulation box is 153. 6cr x 79.06cr x 79.06a and the chain 
length N is 96. It yields the monomer concentration C m equal to 0.0001cr~ 3 , which is about 
12 mM in real unit. This C m describes a dilute PE solution. The salt concentration C s is 
set to the equivalence point C* = C m /3. It has been shown that at C s = C*, the trivalent 
counterions form complexes with the chain, lar gely neutralizing the bare chain charge; the 



chain collapses into a globule structure 26j, [271, |35|, |36[. $ and \E\ are varied over a wide 



range of value to study their effects on the properties of the PE system. Langevin equation 
is integrated by Verlet algorithm with a time step At equal to 0.005r j^]. A pre-run phase 
takes about 10 7 time steps to bring the system into a stationary state for each |J5|)-case, 
followed by a production-run phase of 10 8 time steps to cumulate data for analysis. The 
hydrodynamic interaction is not included in this work because it is largely screened out 



under typical electrophoretic conditions [19 



38 



point 



. Recent study clearly demonstrated this 
40(. To simplify the notation, physical quantities in the following text will be reported 
in (a, m, k B T, e)-unit system if not specially mentioned. For e.g., the strength of electric 
field will be described in unit fcsT/(ecr). 
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III. RESULTS AND DISCUSSIONS 



A. Degree of chain unfolding 



We first revisited a reference case in which the applied field is a DC electric field. It is 
the limiting case of an AC field with infinite \F The degree of chain unfolding D s , defined 
as the ratio of the end-to-end distance R e of chain to the contour length L c , is shown in 
Fig. [T](a) against the DC field strength E dc . 
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FIG. 1: (a) D s = R e /L c in DC fields as a function of E dc ; (b) D s as a function of ^ (in unit At) 
in AC fields of different \E\. The value of \E\ is indicated in the inset of the figure. The symbol 
'X' denotes the inflection point on each corresponding curve. The error bar of data, if not shown, 
is smaller than the symbol size of data in this paper. 



We can see that D s shows a sharp transition, occurred at E dc = 0.22. Below E dc , the 
PE chain stays in a collapsed state, in which condensed trivalent counterions tightly bind 
monomers together, resulting in a chain size similar to the unperturbed size in zero field. The 
structure of the collapsed chain is a spherical globule. When E dc > E dc , the chain unfolds 
abruptly. The DC field is strong enough to drive the condensed trivalent counterions towards 
one side of the chain globule and even more, strips them off the chain. The chain is hence 
elongated and aligned in the field direction in the process. The degree of chain unfolding 
can be as large as 90% of the contour length when E dc is large enough. The detailed study 
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including the effect of salt concentration and valence on chain unfolding as well as the scaling 



law for E^c against N in DC fields can be found in Refs. 3l|, |32| . An independent simulation 
done recently by other group 41[ also confirms our results. 

We then studied the case with AC fields. The degree of chain unfolding is now plotted 
against the field period \1/ in Fig. 1(b) for different field strength \E\. The results show that 
D s takes the value of zero field, 0.07, over the whole range of \I/ when \E\ < E^ c . In this 
situation, the field strength is too weak to overcome the binding of the condensed trivalent 
counterions to unfold the chain. When \E\ > E^ c , chain unfolding can take place and D s is 
a sigmoidal function of When \1/ is small, the applied field alters very frequently. The 
duration of a half period of the field, within which the square-wave field acts as a DC one 
pointing to either +x or — x direction, is not long enough either for the chain to complete 
its conformational change or for the condensed counterions to migrate to stable positions. 
The degree of chain unfolding is thus small. When \l/ is sufficiently large, the deformation 
and the polarization of the PE complex can be completed in a half period. The chain is 
unfolded and the maximal degree of chain unfolding takes a value corresponding to the 
limiting value in the referenced DC field. In the figure, the limiting value has been indicated 
in dashed line. The reduction of the degree of chain unfolding in the region of small \I/ is 



in agreement with experiments 23|, |24J; recent simulation 41] also reported a consistent 



behavior. Moreover, we observed that the transition of D s is not really sharp but occurs 
over a range of period In order to characterize this behavior quantitatively, we define 
the transition point at \1/ = \l/ c , which is the inflection point of the D s curve. The inflection 
point has been indicated by a symbol 'X' on each curve in Fig. H(b). We saw immediately 
that \l/ c depends strongly on \E\ and the value falls in a range from 3 x 10 3 At to 7 x 10 5 At, 
or equivalently, from 22.5 ps to 5.25 ns. Therefore, the collapsed chain is unfolded when the 
AC field frequency is smaller than the order of tera- or giga-hertz. Detailed study showing 
that \l/ c is directly related to the characteristic time to fully polarize the PE chain will be 
presented later. Since chain unfolding occurs in a small frequency region, experiments can 
verify this phenomenon easily, provided that the AC field strength is stronger than E^ c . As 
having been shown in our previous DC simulations {31! 32]. E^ c can be largely reduced to 
a value of few hundreds volt per cm if a long PE chain of order, for e.g., N = 10 6 , is used. 
This strength of AC field is accessible by experiments. 
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B. Simulation Snapshots 



In order to have the picture of chain conformation in mind, we show in Fig.[2]the snapshots 
of PE chain in AC fields with the field strength equal to 0.3. 




The PE chain is represented by the yellow-colored bead-spring chain. The trivalent counterions, 
the monovalent counterions, and the coions are represented by the red, green and white beads, 
respectively. The change of the field direction is indicated on the left side of the figure as a 
function of field cycle t/^f and the evolution of the chain conformation can be seen in each panel, 
following the t/^-axis from bottom to top. 

Three field periods are presented in the panels: (a) $ = 10 3 At, (b) f = 8 x 10 4 At, and 
(c) \I/ = 10 6 At, which are, respectively, smaller than, about equal to, and larger than the 
critical period \fr c (refer to Fig. [T]). In each panel, the evolution of chain shape against a 
field cycle is shown, following the t/^/-axis. We can see that for the case of small period 
(Fig. EJ^a)), the chain is not unfolded in the AC field. It exhibits a globule structure. This 
is because the acting time in a half-cycle of the AC field is too short for the chain to 
response. On the other hand, for the case of large period as shown in Fig. EJ^c), the acting 
time is long enough. The chain is stretched and displays an elongated structure aligned 
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parallel to the field direction. If we look carefully into the snapshots, we can find that the 
trivalent counterions are not uniformly condensed on the chain backbone in this case. They 
follow the field direction and accumulates more on one of the chain end. This non-uniform 
distribution results in a chain polarization. When the AC field switches its direction, the 
stability of the system is suddenly destroyed and the condensed trivalent ions migrate toward 
the other chain end. The chain partially shrinks and then regains its size soon after the new 
stability is established. For the case with the field period close to \l/ c (Panel (b) of Fig. [2]), 
an astonishing, tadpole-like structure of chain is observed, just before switching of the AC 
field direction. The head of the tadpole structure is derived from the accumulated trivalent 
counterions binding together with the chain near an end. And the tail is formed by the 
other side of the chain with less counterions condensed on it. Because the chain is partially 
unfolded, the degree of polarization is smaller than the one in Panel (c). The change of 
the AC field direction drives the trivalent counterions moving away from the tadpole head 
toward the tail. As a passage, the head shrinks and the tail grows. Eventually, the new 
accumulation of the trivalent counterions creates a new head on the other side of the chain. 
The new tadpole structure thus inverts its heading direction and waits for the next field 
alternating. 



C. Theory of chain unfolding 

The dependence of \l/ c on \E\ can be captured by the following charge relaxation theory. 
We assume that a PE chain is unfolded in an AC field if the polarization energy of the chain 
complex, W ac = a(u)\E\ 2 /2, is larger than some critical energy, where u is the field frequency 
and a{u) is the polarizability. This critical energy can be estimated by W£ c = a(0)E*£/2, 
the polarization energy to unfold the chain in the critical DC field E^ c . It is known that for 
a conducting Maxwell- Wagner sphere of arbitrary radius R in an AC field, the polarizability 
can be expressed by a(u) = 4ire m K(u)R 3 where 



K{u) 



e p {u) - e m {u) 



e p (uj) + 2e m {uj) 

is the Clausius-Mossotti factor, and e m and e p are the cu-dependent complex permittivity of 



the solvent medium and the PE complex, respectively 42|. This Clausius-Mossotti factor 



arises due to both the dielectric polarization, resulting from the permittivity difference across 
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the sphere-medium interface, and the space charge accumulation at the interface, derived 
from the conductivity gradient and the corresponding ion flux discontinuity at that location. 
While counterion migration occurs across the entire collapsed chain, charge polarization 
occurs only at the two poles. It involves the dissociation and condensation of the migrating 
counterions presented within a Debye-layer neighborhood of the interface. As a result, the 
dipole formation time is the usual charge relaxation time at the interface. When the resulting 
dipole moment aligns parallel to the electric field, a net stretching force is imposed on the 
chain, corresponding to the polarization energy. The chain is unfolded. The condensed ions 
hence play a role to the PE chain, similar to the one of the bound surface electrons to a 
dielectric object. Therefore, the polarization dynamics of the system can be studied by the 



method of electrical impedance 



42 



42| 



We model the dielectric response for the PE complex by a resistor-capacitor (RC) circuit 
in series. The complex permittivity reads as 



e p {u) - 



1 + iuJSp/ a p 1 + iujT rC)P 
The response function for the solvent is modeled by a RC circuit in parallel and the complex 
permittivity is 



Here e p and a p (or e m and a m ) are, respectively, the permittivity and the ion conductivity in 
the PE complex (or in the solvent), and r rCiP = e p /a p and r rCjm = e m /a m . Since the trivalent 
counterions condensed quite tightly on the chain, the migration of these ions from one side 
of the PE complex to the other side in an electric field is similar to the electron migration 
from one plate of a capacitor to the other plate in a voltage source. The ion migration also 
suffers a resistance force against the applied electric field, a situation which is similar to 
put a resistor in the circuit. Therefore, we model the PE chain phase by an RC circuit in 
series. Because the length scale across the phase in the field direction corresponds to both 
the capacitor gap and the resistor length, the ratio r rcp is the characteristic time to move 
the counterions from one side of the PE chain to the other side. It corresponds to the RC 
time of the circuit. In the solvent phase, the ions travel through the whole system. Since the 
traveling speed is not fast, the medium is a poor conductor for ions. Therefore, we regard 
the solvent as a leaky dielectric, modeled by an RC circuit in parallel. A similar argument 
shows that the ratio r rCjm is the RC time for the solvent. 
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Due to the linear dependence of the conductivity on the ionic strength, the RC time for 
solvent is also the charge relaxation time of electrolyte, which is equal to the squared Debye 
screening length divided by the diffusivity of the dominant ions. For the PE phase, this 
connection is typically not true because the ions conduct through nanoporous structures 
formed by the twisted chain in the phase, rather than simply through a bulk solution. 
The RC time for the PE phase is important in determination of the critical frequency. 
At high frequency, the chain polarization is due to dielectric polarization and the value is 
small. At low frequency, the "capacitor" is fully charged through conductive polarization 
by accumulating counterions at one pole of the PE globule, producing a large polarization. 
This polarization can grow even further via chain unfolding if the applied electric field is 
strong. 

We assume that the RC time for PE is a constant at the critical points to unfold the 
chain. Based upon the above assumption for the polarization energy, the critical frequency 
uj c to unfold a chain is determined by the criterion (\E\/E^ C ) 2 = Re[K(0)]/Re[K(u c )}. Let 
us consider the case of e p ~ e m . We obtain 

r-n _ (Q T rc,p T rc,m + ^r^rni^c ^ 

4 + (4r 2 cp — 2T rCi pT rCtm )UJ 2 + 4n~r c,p T r 0,171^ c 

where H = (|.E|/i^ c ) 2 — 1 is the excess of the squared electric field. When u c is small, the 
equation gives 5 ~ ( 

+ j T rc,m) u c + C( w c) ; which is a linear function of The 
figure H vs. w 2 is plotted in Fig. [3] where uj c is calculated by 2ti/^ c . 

The result does show a linear dependence of 5 in the region of small a; 2 , which follows the 
prediction of our theory. Least square fit by Eq. (JH]) using Levenberg-Marquardt algorithm 
yields T rCjP = 0.393(8) and T rc>m = 21.2(2). In our simulations, the length unit a and 
time unit r are 2.4A and 1.5ps, respectively, while mapped to a real system. Since the 
trivalent cations condense on the chain and largely neutralize the chain charge, most of the 
monovalent cations and anions are presented in the bulk solution. Hence, the concentrations 
of the cations and the anions in the solution are both about 0.0001 (cr~ 3 ), or equivalently, 
12 niM. The permittivity of water in our system is 

e 2 _L_ e 2 

r AnksTXB 127r k^To 

Using this value of e e r (1/127T in our unit system) for s m , the theoretical interpretation of our 
coarse-grained simulation of chain stretching yields the solvent conductivity a m (= £ m /r rCjm ) 
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FIG. 3: E 

curve is the result of fitting by Eq. ([5]). 



equal to 3.5 x 10 , or equivalent ly, 0.59 S 
with the one of sodium chloride in water 



in 



44J. This solvent conductivity agrees fairly 



451 ] . providing a self-consistent support for our 



theory. Here we have corrected the conductivity by dividing the value by a factor 35. It 
is because the friction coefficient in Eq. (jlj) has been set to a small value Q = l.Or -1 in 
our study, which is about 35 times smaller than in water, for the purpose to accelerate 
the simulations. Therefore, to map diffusion-related quantities such as conductivity and 
relaxation time, we have to correct the obtained values by the factor 35. 

The typical RC relaxation time in exp eriments ranges from milliseconds to nanoseconds, 
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depending on the ionic strength 
I s = 0.5 ((+1) 2 x 10- 4 + (-1) 2 x 10" 4 



42]. The ionic strength in our bulk solution is about 
10 -4 (er~ 3 ), which gives a Debye screening length 



equal to In = (47rA# • 2/ s )~ 1 / 2 ~ 2.76 nm. Using the typical value of ion diffusivity in water, 
Di ~ 10~ 9 m 2 / s [46j , we expect the relaxation time in water with a value of about ij^/Di ~ 7.6 
ns. Inside the PE globule, the ions diffuse through the nanoporous structure formed by the 
twisted chain. The size of the nanopores is about a, the diameter of the ions. Therefore, 
the relaxation time in PE phase is estimated by u 2 / Di and the value is 58 ps. The RC time, 
r rC)V and x rc m obtained from Fig. [3l is 20.6 ps and 1.1 ns, respectively, after multiplying the 
correction factor. These values are in fair accordance with the theoretical estimation. We 
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remark that the RC time of highly-concentrated nanocolloid or molecule suspensions can be 
measured from dielectric relaxation spectroscopy or impedance spectroscopy, and the value 
for a single nanocolloid or molecule is inferred by an averaging theory 42J. Recently, cross- 
over frequencies of dielectrophoresis have been used to estimate the RC time of individual 
blood cells and nanocolloids j^, (47)]. A very recent study has demonstrated an approach to 
determine the RC time of short PE using fluorescent correlation spectroscopy, based upon 
the coil-globule transition in AC fields [251 ] . 



D. Time variation of dipole moment and chain unfolding 

The critical frequency u c obtained in Sec. III.C is independent of the radius of gyration 
of the chain, as the Maxwell- Wagner relaxation time describes the relaxation/dissociation 
of ions within a Debye layer from the interface. This physical picture breakdowns when the 
chain begins to uncoil, such that the trivalent counterion no longer bridges two segments 
of the polyelectrolyte. Internal polarization could also develop at a different time scale. In 
order to understand the coupling between chain unfolding dynamics and interfacial charge 
relaxation/dissociation dynamics, we study here how p x and D s evolves in cycles (t/^f) of 
the AC field for different values of p x is the induced dipole moment p of the PE complex 
in the field direction. The PE complex, by definition, consists of the chain itself and the 
ions condensed inside a tube region of radius r t = Ab surrounding the chain backbone, p is 
calculated by J2i Zi^ifi — r C m) where i runs over all the particles in the PE complex and r* cm 
is the center of mass of the chain. The results for \E\ = 0.3 are presented in Fig. HI 

We saw in Fig. H)^a) that p x grows exponentially with time in every half cycle of the 
square-wave field and switches its sign simultaneously with the field. The growth of p x 
saturates to some value if \I> is large enough, for example, ^ = 10 6 At. It suggests that 
there is a relaxation time to polarize the PE complex in an electric field. Fig. SJb) shows 
that for large the chain shrinks right after the field switches the direction, but the chain 
size does not reduce down to the zero-field value. It regains its size in the corresponding DC 
field after p x reaches the saturation value. On the other hand, for small the variation 
of D s is not significant. To go further, we drew the parametric plot of p x and D s with 
time in Fig. SJc). We observed that the trajectory points travel along a butterfly curve, 
from the tip of one wing to the tip of the other wing, after alternating the field direction, 
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FIG. 4: Time variation of (a) p x and (b) against the cycle of AC field, (c) Parametric plot of 
p x and D s with time. The field strength is \E\ = 0.3. The value of \I/ (in unit At) is indicated near 
the associated curve and data. 



and then wander near the region of wing tip waiting for the next switch of field direction, 
as sketched in the inset of the figure. The wing bottom occurs at non-zero p x , showing 
that the variation of D s falls behind the p x . Therefore, the chain elongation is slaved to 
the fast dynamics of the interfacial trivalent counterion dissociation/relaxation time which 
dominates the polarization. Moreover, a beautiful shuttlecock pattern is formed when the 
ensemble of the data at different field frequencies is looked together in the p x -D s plot. 



E. Polarization time and fluctuation time 

To quantify the findings with our theory, we calculated the polarization time by fitting 
exponentially the growth of p x in a half cycle of field with large The mean polarization 
time t p versus the critical period \I/ C to unfold the chain is shown in Fig|5^a). 

We found that r p depends linearly on with slope equal to 0.27(1). This result inter- 
estingly connects the critical points with the polarization in the supercritical region where 
the chain is fully unfolded. The physics can be understood by the following simple picture. 
Because the chain structure is linear in the supercritical region (\l/ > \l/ c ), the polarization 
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FIG. 5: (a) polarization time r p vs. critical period Vl/ C , and (b) t p vs. field strength \E\. In panel 
(b), fluctuation time Tf is also plotted as a function of E^c at the same scale. The values of t p , Tf 
and \& c are reported in unit At. 



time t p can be estimated by the relaxation time R/vd where R is the relaxation distance 
of trivalent counterions migrating through the chain to establish a new polarization when 
the applied field switches the direction and is the drifting velocity equal to n e \E\. In our 
simulations, /i e ~ 3 for the trivalent counterions and \E\ ~ 32E^ c cu c for small u c according 
to Eq. ([5]). The slope 0.27 yields a relaxation distance of R ~ 36 (about 30% of the chain 
contour length), which is a reasonable value in view of the profile of charge distribution along 
the chain. Fig[5^b) shows r p as a function of the AC field strength \E\. One can see that the 
stronger the field strength, the shorter the polarization time. It is because ions and PEs are 
dragged at faster speeds when the applied field is stronger. Moreover, detailed study showed 
that in the region of strong field \E\ 3> E* dc , t p is inversely proportional to \E\. Also, in the 
limit where \E\ tends toward E^, the product of and IE"! 2 — E^ is a constant. The two 
behaviors can be deduced easily from our theory, known that \I/ C and t p are identical within 
a factor. Our simulations are in good agreement with experiments of field-induced decon- 



densation 



2l|, in which the reciprocal of the relaxation time of decondensation exhibits a 



linear relation to the field strength when the field is strong, and deviates from the relation 
when the field is weak. 
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Recently, simulations done by Liu et al. 4l| suggested that the intrinsic chain relaxation 
frequency is the upper bound of the critical AC frequency to stretch a PE chain. This bound 
is quite crude in light of the current result. To further verify the importance of the dissocia- 
tion/condensation of the interfacial trivalent counterions, we calculated the autocorrelation 
function of D s , defined by 

{DD - (D s ) 2 ' 

in different DC electric fields. The correlation time 77, calculated by fitting R(t) function 
exponentially, represents physically the fluctuation time of chain size under the action of a 
DC electric field. The results have been plotted in Fig 0(b) as a function of the field strength 
E,i c at the same scale. We found that 77 is a constant when E& c < E dc . A large jump occurs 
at Edc = E* dc and then, 77 decreases inversely with increasing E& c . A striking finding is that 
77 is basically coincident with t p in the region E^ > E* dc . We know that the PE chain is 
in a dynamic equilibrium with the ions surrounding it. The dissociation or condensation 
of a counterion away or onto the chain will locally break the chain stability and cause the 
fluctuation of the chain size. The coincidence between these characteristic times shows the 
importance of the ion fluctuations on chain unfolding. The critical frequency to stretch a 
chain in a square- wave AC field is mainly determined by the dissociation/relaxation of the 
condensed interfacial trivalent counterions, which takes the value of the fluctuation time of 
chain size in the counterpart DC electric field. 



IV. SUMMARY 



We have revealed the key conditions to unfold collapsed PE chains in square-wave AC 
electric fields: the field strength must be larger than the critical DC field to stretch the 
chain and the frequency is smaller than the inverse of the counterion dissociation/relaxation 
time under a quasi-DC field. Moreover, when both conditions are satisfied, there is a strong 
correlation between the critical field strength and the field frequency. A simple scaling theory, 
based upon Maxwell- Wagner dielectric theory, can quantitatively capture this critical field- 
frequency correlation for PE unfolding in AC fields. 

This material is based upon work supported by the National Science Council, the Re- 
public of China, under Grant No. NSC 97-2112-M-007-007-MY3. Computing resources are 
supported by the National Center for High-performance Computing. HCC is supported by 
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